NONLINEAR SCHRODINGER EQUATIONS WITH 
MULTIPLE- WELL POTENTIAL 

ANDREA SACCHETTI 

O • 

^». 1 ^ Abstract. We consider the stationary solutions for a class of Schrodinger 

equations with a A^-well potential and a nonlinear perturbation. By means 
of semiclassical techniques we prove that the dominant term of the ground 
state solutions is described by a Af-dimensional Hamiltonian system, where 
the coupling term among the coordinates is a tridiagonal Toeplitz matrix. In 

>J ' , particular we consider the case of Af = 4 wells, where we show the occurrence 

of spontaneous symmetry-breaking bifurcation effect. In particular, in the 
limit of large focusing nonlinearity we prove that the ground state stationary 
solutions consist of TV wavefunctions localized on a single well. 
PACS number(s): 05.45.-a, 02.30.Oz, 03.65.Sq, 03.75.Lm 

Keywords: Nonlinear dynamics, Bifurcation, Semiclassical limit, Bose-Einstein 
condensates in lattices 
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1. Introduction 



In these last years there has been an increasing interest in the study of nonhnear 
Schrodinger equation with an external potential. In fact, many interesting and 
current physical problems may be described by means of such a model, e.g. Bose- 
O^ . Einstein condensates [13) . non-linear optics [6^, semiconductors 12 , and quantum 

l/^ ' chemistry [3j [7], just to mention the most relevant. In particular, the mathematical 

r«^ , research recently focused on the nonlinear Schrodinger equation (hereafter NLS) 

f^ ' with double well potential. One of the most interesting feature of such a model is 

the spontaneous symmetry breaking phenomenon [5j I10[ I14j , and recently a general 
rule in order to classify the kind of bifurcation has been obtained [21 [15] (see also 
[9]). Much less is known for NLS with multiple well potential. So far, few models 
with multiple wells have been considered, e.g. the moden with three wells on a 
regular lattice 8] (where a lattice means a sequence of points displaced along a 
C^ ' straight line) , and the model with four wells on the vertex of a regular square [16) . 

The generic case with N wells has not been yet studied. In fact, multiple-well 
potential represent the effect of small lattice on, e.g., Bose-Einstein condensates; 
furthermore, they are also interesting in order to understand the transition to a 
lattice with infinitely many wells. 

In this paper we consider a NLS with N wells displaced on a regular lattice (even 
if the present analysis can be easily extended to the case on wells displaced on a 
regular grid, as discussed in an explicit example in Remark |4l). We'll show that 
such a problem can be reduced, up to a remainder term, to a finite dimensional 
system; to this end, instead of using some kind of Galerkin decomposition as in [8], 
we assume to be in the sem,iclassical limit. In such a way we can make use of some 
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2 ANDREA SACCHETTI 

powerful! results of the semiclassical analysis [4^ . By means of such results and by 
making us of the Lyapunov-Schmidt reduction scheme we estimate the remainder 
term. 

The finite-dimensional system we obtain is almost decoupled, in the sense that 
the coupling term which represents the interaction among the adjacent wells is 
associated to a tridiagonal Toeplitz matrix. Furthermore, it can be written in 
Hamiltonian form, where one of the coordinate is a cyclic coordinate. 

We then consider in details the case of TV = 4 wells and we study the bifurcation 
picture when the strength of the nonlinear perturbation increases. As in the 
double well model we can see that the ground state stationary symmetric solution 
bifurcates giving arise to 4 stationary solutions fully localized on a single well, and 
the kind of bifurcation satisfies the same rule as in the double well model. Actually, 
such a result may be generalized to any number N of wells by means of a simple 
asymptotics argument. 

2. Description of the model 
Here, we consider the nonlinear Schrodinger (hereafter NLS) equations 

^h^^HoyJ + eg{xM'-i:, ||V(-,i)|| = 1, (1) 

where e G R and || • || denotes the L"^ {M.'^ , dx) norm, 

Ho = -£a + K A^t^ (2) 

is the linear Hamiltonian with a multiple-well potential V{x), and g{x)\il;\'^'^ is a 
nonlinear perturbation. For the sake of definiteness we assume the units such that 
2m = 1. The semiclassical parameter h is such that h <^ 1. 

Here, we introduce the assumptions on the multiple-well potential V and we 
collect some semiclassical results on the linear operator Hq. 

Hypothesis 1. Let v{x) € C^{M.'^) be a smooth com,pact support function with a 
non degenerate minimum value at x = 0: 

v{x) > v,mn = v(0), Va; e R^ x^ 0. (3) 

We consider multiple-well potentials of the kind 

N 

V{x)^Y.'^{x-x,) (4) 



for some N > 1, where 



j=i 



^^,0,. 



and £ > 2r, where r > is such that 

C C [-r, +r] X R'^^i 

and where C is the compact support of v{x). 

Hence, the multi-well potential V{x) has exactly N non degenerate minima at 
x^Xj, j ^1,2,...,N. 



NON LINEAR MULTIPLE-WELL 3 

Remark 1. It is a well known fact [2 that the Cauchy problem (CJ) is globally 
well-posedness for any e S R small enough provided that 

( +CXD if d <2 
In such a case the conservation of the norm of^p(x,t) and of the energy 

(T + 1 

follows; furthermore we also have a priori estimate 

m-.tf^\\<ch-'^ (5) 

for some positive constants C and m. 

3. Analysis of the linear Schrodinger equation 

Now, making use of semiclassical analysis [4 we look for the ground state of the 
linear Schrodinger equation 

Let dA{x,y) be the Agmon distance between two points x and y, let 

Sq = midA{xj,Xi); 

then, by construction of the potential V{x), it turns out that 

So ^dA{xj,Xj+i), j ^1,2,...,N -1, (6) 

and 

So <dA{xj,Xi) if |i-j|>l. (7) 

Now, let Hd be the Dirichlet realization of 

Hd = -h^A + V (8) 

on the ball Bs{0) with center at a; = and radius S > 2So- Since the bottom 
of v{x) is not degenerate, then the Dirichlet problem associated to the single-well 
trapping potential v{x) has spectrum a{IIr)) with ground state 

d 

XD = Y.^h + oih^), 

where 2/Xj are the positive eigenvalues of the Hessian matrix v"{0), such that 

dist [XD,cr{HD) \ {Xd}] > 2Ch 

for some C > 0; the associated normalized eigenvector ^^{x) is localized in a 
neighborhood of a: = and it exponentially decreases as O (^fi-"^g-dAix)/h'^ -fQj. 
some TO > and where dA{x) is the Agmon distance between x and the point 
x = 0. 

The spectrum cr{Ho) of Ho contains exactly N eigenvalues Aj, j = 1, 2, . . . , iV, 
such that 

A, -Xd^ Oie^p"") 

for any Q < p < So] furthermore 

dist [{A,}f^i,a(ifo) \ {A,}f^i] > Ch. 
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Let F be the eigenspace spanned by the eigenvectors ipj associated to the eigenvalues 
\j. Then, the restriction Ho\f of Hq to the subspace F can be represented in the 
basis of orthonormahzed vectors 0, such that 



0,(x) - ^,{x) - O (e-"/'') , j = 1, 2, . . . , TV, ifi.ix) = ^oix - x,), 



(9) 



for any fixed < p < 6*0; the eigenvector (j)j is locahzed in a neighborhood of the 
minima points Xj. More precisely, Hq\f in the basis (f>j{x), j = 1,2,...,A^, is 
represented by the N x N matrix (see, e.g. Theorem 4.3.4 by [4]) 



A,3lw + K,)f,=i+0(?i°°)e-^«/'' 



where 



and where 



Wi, 



0, if N-jIt^i, 



(10) 
(11) 



is independent of i = 1, . . . , iV and it is such that (see Theorem 4.4.6 by [4]) 



(12) 



Remark 2. Let i — \, . . . ,N — \ he fixed and let Vt he an open set with smooth 
houndary such that Xi € Q, and Xj ^ J7 for i ^ j , let 

Ef ^ {x : dA{xi,x)+dA{xi+i^x) < SQ + a} 



for some a > sufficiently small. Let T i = dfl D E.^ 



(a) 



the 



dipD(x-x^+i) dt/joix-x^ 

1pD{X - Xi) ^d{X - Xi+i)- 



dn 



dn 



dSr- (13) 



The dominant term of j3 is independent of a. 

In particular, in dimension one, i.e. d = 1, then it turns out that 

Collecting all these results then we can conclude that 

Lemma 1. Let F be the eigenspace spanned hy the eigenvectors 4'j associated to 
the eigenvalues Xj of Hq. Then, the restriction Ho\f of Hq to the suhspace F can 
he represented in the basis of vectors 4>j, localized on the j—th well and satisfying 
(0), by the N x N tridiagonal Toeplitz matrix 



ca\ -So/h 



where 



T + 0{h°°)e 



T^-/3T+Xd1 



N 



(14) 
(15) 



NON LINEAR MULTIPLE-WELL 



r = 



/ 


1 
















\ 


1 





1 



















1 





1 



















1 





1 





















1 





1 



















1 





1 


\o 
















1 


0/ 



(16) 



and satisfying /11S\) for some 



where /3 is the positive real number given by U 
C>0. 

From (J14p it turns out that the eigenvalues Xj are given, up to a sniaU correction 
0{h°°)e~^°'^, by the eigenvalues fij of the matrix T. To this end we recall that 
the A^ eigenvalues of the tridiagonal Toeplitz matrix T are given by |llj 



fij ^ Xd — "2/3 cos 



N + 1 



j = l,2,...,iV, 



with associated eigenvectors 



(«.) 



3)k 



kJTT 



In order to normalize the eigenvector we remark that 

kJTT \ N +1 



N 



Y^sii 



fc=i 



^N+lJ 2 

From this fact and from Lemma [1] we then conclude that 
Lemma 2. The first N eigenvalues Xj of Hq are given by 



Xj = A_D — 2/3 cos 



0{h^)e- 



J = l,2,. 



and the associated normalized eigenvector are given by 

N 

i^jix) = J2 c^j,kM^) + 0(ft°°)e-^»/'' , J = 1, 2, . . . , TV, 



,N, 



(17) 



fc=i 



where 



ctj,k = afcj = 



2 / kj-n: 

■ sin f 



iV+1 



A^ + l 



and ^kix) ^ ^d{x ~ Xk). 



Remark 3. In the case of two wells, i.e. N — 2, then the two eigenvalues are 
given by 

Xi = Xd — (3 and X2 = Xd + P 

and the matrix A = {oljiS) has the form (in agreement with 15J 

V2/I 1 
^~ 2 U -1 
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Figure 1. Double- well potential (left) and the two symmetric 
and antisymmetrical eigenvectors (right). 



In the case of three wells, i.e. N = 3, then 

Ai = Au — v2/3 , A2 — \d and 
and (in agreement with [S]^ 

1 %/2 



A3 — Ar 



V2I3 




In the case of four wells, i.e. N = 4, then 



Ai 



2/3 COS (- 



A2=A 



D 



,'27r 
2^ cos ( — 



\d + 2/3 cos 



2tt 



A4 = Xd + 2/3 cos f-] 



and 



A 



I sin(7r/5) 
sin(27r/5) 
sin(27r/5) 



sin(27r/5) 

sin(7r/5) 

— sin(7r/5) 



sin(27r/5) 

— sin(7r/5) 

— sin(7r/5) 



sin(7r/5) \ 
- sin(27r/5) 
sin(27r/5) 



\ sin(7r/5) - sin(27r/5) sin(27r/5) - sin(7r/5) / 



For instance, see Figures 1, 2 and 3 for the one- dimensional N-wells problem 
with, respectively, N — 2, N — 3 and N = 4; where the wavefunctions ipoix) and 
ipj (x) can be chosen to be real-valued functions. 

Remark 4. We can immediately extend such an analysis to multiple-well potentials 
of the form V{x) = ^ v{x — Xj) where Xj e W^ are points on a regular grid. We 
study, for argument sake's, the model in M^ considered by [Id] where the potential 
V{x) has 4 wells with minima on the points 

XI = (1, 1); X2 = (1, -1); 0:3 - (-1, 1); .T4 = (-1, -1) ■ 
In such a case we have that (recalling that Wi_j — Wj^i) 

W'1,2 = Wi^3 = W2,4 = W4,3 = ~/3 and Wi^i = 0, Wi,4 = W2,3 = 
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Figure 2. Three- well potential (left) and the three symmetric 
and antisymmetrical eigenvectors (right). 




'V \f 



^r 



JV_ 



Figure 3. Four-well potential (left) and the four symmetric and 
antisymmetrical eigenvectors (right). 



and the matrix T takes the form 

/ Xd -(3 ^13 \ 

^^ -P Xd -/3 

-(3 Xd -/S 

V -13 ^13 Xd J 

with eigenvalues (and associated normalized eigenvectors {vj)k) 

•i 1 1 r 

2' 2' 2' 2. 



Ml 

M2 
M3 

Ml 



Ad -2/3, vi 
Xd , V2 = 
An , W3 = ( 0, 
Az) - 2/3 , wi = 



71-°'°-71 



1 1 



,0 



V2'V2 
1111 

2'"2'^2' 2 
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in agreement with |16]. 

4. The A^-mode approximation for the NLS equation 
Let ip be the normalized solution of the NLS equation ^ written in the form 

N 

By substituting (fT8)) into ([T]) and, projecting on the eigenspaces spanned by the 
eigenvectors tpij and on the orthogonal eigenspace, we obtain the following system 
of differential equations: 

( ihcj = A,c,+e(^,(-),5(-)lV'(-,0l''V'(-,<)) .,q^ 

\ thijjc = i?oV'c + ene(7(x)|V'(x,i)|2-v(x,t) ^ ' 

where we set 

N 

n = ^(-0j, •)V'j and nc = i-n. 

By substituting (|T7| in (|T9| we obtain 

fe 

fc 

= (Az3 + ojj)cj + e^ aJ''^'Cfc|dfcp'"4 + r^ 
fc 
where we set 

N 

dfc:=dfe(i)=^c,(i)a,-fc, (20) 

i=i 

^.--2/3cos(^^) and C, = {^l+\-),g{-)^l+\-)) , 

since ((p^, I'^mp'^'yam) = 0{h°°)e~^°^^ for m ^ fc and a^^fc = aj^k = o^kj, and where 
we make use of the a priori estimate ([5]) of the norm of |'0P'^- We should underline 
that, by construction and by Lemma [2l it follows that 

N 

fe=i 

If we denote by 

A = {aj^k), C = (Cfe), d = (rffe), A = diag(Aj), r = {rk)k , 
then the above equation takes the form (with abuse of notation) 

ihA-^d = AA-^d + eCA-^ldf^d + r 
that is 

ihdk =- {Td)k +eCk\dk\^^dk+r (21) 
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since 

T = AKA-^ 

and where we set 

C = ACA-^ and f ^ Ar. 

Definition 1. We call A^-mode approximation for the NLS equation the system 
of ODEs obtained by neglecting the remainder term r 

ihdk = {Td)k + eCk\dkf"dk, k^l,2,...,N, (22) 

where Ck are real-valued constant, and with the normalization condition 

N 

J2\dk{t)\' = l. (23) 

fe=i 

Remark 5. It has been proved the validity of the N-mode approximation for times 
of the order /3~^ by pj. Concerning the stationary solutions has been proved by 
[3j that the 2-mode approximation gives the stationary solutions for the NLS, up 
to an exponentially small error, furthermore the orbital stability of the stationary 
solutions is discussed; the same argument may apply to the N-mode approximation 
for any N > 2 proving that the stationary solutions of equations i2S\) and \2S\) give, 
up to an exponentially small error 0{e^P'^), for any < p < Sq, the stationary 
solution of the NLS {IJ). However, we don't dwell here on these details. 

Remark 6. For instance, in the case of two symmetric wells, i.e., N = 2 then ^27 
takes the form (in agreement with [15j J 

ihdi = Audi — /3(i2 + eCildipdi 



ihd2 = \Dd2 - Ml + eC2\d2\^d 



2 



In the case of three symmetric wells, i.e., N — 3 then \22(l takes the form (in 
agreement with [8]j 

ihdi — Xodi— I3d2 + eCi\di\'^di 

ihd2 = \Dd2- Pdi- l3d3 + eC2\d2\'^d2 



ihds = Xods- (3d2 + eC3\d3\^d 



3 



Finally, in the case of four symmetric wells, i.e., N ~ A then \2'/3jl takes the form 
ihdi = Xodi — I3d2 + eCi\di\'^di 



(24) 



ihd2 = \Dd2- I3di- pd3 + eC2\d2\^d2 
ihds = Xnds- (3d2- Pdi + eCsldsl^ds 
ihd/i = \odi — I3d^ + €Ci\di\'^d/i 

4.1. Hamiltonian form of the A^-mode approximation. If we set 

4 = V*e*'N Ik^ [0,1], Ok^ [0,2^), 

then, by means of a straightforward calculation, it follows that (|22]) takes the Hamil- 
tonian form 



^9fe = Hf = -~'^PY.j=i Tj,ky^WqjSm{0j - 



huk — — a — — — 

" aqk 



^D-/3 E =1 Tj,kJ^ cos(0, - dk) + eCkql 



(25) 
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with Hamiltonian function 

H := 'H{qi,...,qN,0l,---,0N) 

N N 1 ^ 

= Xn^Qk- P ^ Tkj cos{9j - Ok)y/qJq^ + <^-—j Yl ^'^'ik^^ 

k=l j:k=l k=l 

N N~l N 

= Az5 ^ (?fe - 2/? ^ cos(0fe+i - Okj^q^T^ + (—-^ Y^ Ckql^^ . (26) 



k=\ k=l fc=l 

The normahzation condition (j23p takes the form 

N 

E* = l' (27) 

fc=i 

furthermore the Hamihonian function "H is an integral of motion; i.e. 

n[qi{t), . . . ,qN{t),9i{t), . . . ,eN{t)] ^ const.. (28) 

4.2. Reduced Hamiltonian. We make use of the fact that ^j. g^ = 1 in order 
to reduce from A^ to iV — 1 the degree of freedom of the Hamihonian system (|26| . 
We consider the canonical transformation {q,9) — >■ (Q,6) defined as 

h 

Qh = Y qk 

fc=i 

where Qi G [0, 1] and Qh E [0, 1 — Qh-i], h = 2, . . . ,N ~ 1. The inverse transfor- 
mation is defined as 

qi = Qi and qf^:^ Q^- Qh-i , h^2,...,N. 
The associated transformation on the conjugate variable is then given by 

&N = On and 9,, = 9,, - 9h+i , h = 2,...,N, 
with inverse 

N 



E^^ 



k=h 



In the coordinates {Q, 6) the Hamiltonian system takes the form 



flQk — QQ^ 



96^ , fc = l,...,iV, (29) 

dQk 

where the new Hamiltonian denoted by /C is given by 

K = XdQn - 2/3cos(ei)v/(Q2 - Qi)Qi - 2/3 E cos(efc) V(Qfe+i - Qk){Qk - Qk-i) + 



1 

-e- 



>+l 



fc=2 

AT 



fc=2 

It turns out hat 0jv is a cyclic coordinate (indeed Qn = const. — 1), then the 
Hamiltonian system if finally given by 



hQk = f^ 



^ Q^- , k = l,...,N -I, (30) 
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with Hamiltonian function 

/c* = Xd - 2/3cos(ei) V(Q2 - Qi)Oi - 2 J2 cos{ek)ViQk+i - QkKQk - Qk-i) 



Af-2 



fc=2 



-2/3cOs(eA,_i)v/(l - QN-l)iQN-l - Qn-2) + 

1 r w-1 



a+1 



A:=2 



5. Stationary solutions 

Now, we look for the normahzed stationary solutions of the form ip{x,t) = 
e"^*V'(a^)- In. terms of A^-mode approximation ((26|) it consists of looking for the 
solution of the system of equations 



an 

dqk 



(31) 



That is, equation (|5T|) and the normalization condition lead us to the following 
system 



-2sin(6'fc+i - 0k)^/qk+Tqk + '2sm{9k - Ok-i)^/qk-iqk 
-2sin(6'Ar - 6N-i)^qNqN-i 

-Xd - 13 cos{9k+i 
-Xd - (3 cos(6'jv - ^a^-i 



= 
= 
= 



'k)^^+cosi9k-ek-i)^^\+eCkq^, = fuv 



(32) 



-Az5-/3cos(02-^i) 
qi + q2 + ■ ■ ■ + qw 

where A: = 2, . . . ,iV - 2. 



eCigf 



= 1 



Remark 7. We may remark that from the first three equations of f \S^) it follows 
that 

sin(6'fe+i - 6k)^qk+iqk = 0, for any k^l,...,N -1 

Furthermore, from h22fl the stationary solutions where di = and djq = are not 
admitted, hence \dk\ < 1 for any k. Then qi,qN S (0,1) and qk G [0,1) for any 
k — 2, . . . , N — 1. For what concern the phases Ok G [0, 27r) they are defined up to 
a gauge term, then we can always assume that, e.g., 9i — 0. 

5.1. Four wells. The stationary problem in the case of two wells and three wells 
have been already studied [2 [H [15]. We restrict our analysis to the four- well 
case. For the sake of definetness we consider the four-well {N = 4) model where 
we assume that g{x) is a constant function, hence Ck is independent of k and thus 
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we set C := Ck = Ck- In such a case ([32]) takes the form 





sin(02 - 


-Oi)^q2qi 






= 




sin(03 - 


- ^2)1/93 92 






= 




sin(04 - 


~ ^3)1/94(73 






= 




Ad 

/3 


cos(02 - ^1)^ 


/t+Wr 




= ujhl^ 


< 


-ho- 

P 


h ['^"^(^^ - 


-^2)^/9J + cos(6l2 


-0l)^] ^^ql 


= ujh/13 




Ad 


}e [^^"^(^4 


'^3)\/94 + COs(6'3 


-Q2)^^+ml 


= uh/fi 




Ad 

/3 


COS (04 — ^3), 


/f+wf 




^ ujh/p 




^ 91+92 


+ 93 + 94 






= 1 


wher 


e we set 











(33) 



First of all we underline that 92 = (and similarly ^3 = 0) cannot be a solution of 
such a system; indeed if (i2 = is a stationary solution of the 4-mode approximation 
p4l). then ^ reduces to 

iMi — Xddi + eC\di\'^di 

= -/3di - Pds 

ihds = Add3 - Pdi + eCldsl'^ds 

^ ihdi — A(;(i4 — /3(i3 + eC|d4p(i4 

form which follows di = d^ and then ^4 = 0, which is not possible since Remark [T] 

Remark 8. We can extend such an argument to any N : if N is an even integer 
and positive number then the stationary solution are given for qk € (0,1), for any k. 
If N is an odd integer and positive number then the stationary solutions are given 
for gi, qpj € (0, 1) and q^ G [0, 1), fc = 2, . . . , A^ — 1, that is qk — 0, for k ^ 1,N , 
may be admitted values for stationary solutions. 

Since the stationary solutions are given for qk E (0, 1) then it follows that 

ek+i-9k = 0,TT 

and we obtain a family of 8 systems of equations 



-(-ipyf+wf 


= n 






. _-i=[(-l)™^+(-l)^ygj]+^5- = 




, j, m,i = 


1,2 


_ 91 + 92 + 93 + 94 


= n 
= 1 






we set 








^ Xd + hio 









(34) 
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Figure 4. Here we plot the energy il corresponding to sym- 
metrical/antisymmetrical solutions versus rj in the case of cubic 
nonlinearity (i.e. tr = 1). Black line corresponds to the choice 
j — i — 2, and blue line corresponds to the choice j = 1 and ^ = 2, 
that is both graphs correspond to symmetric solutions; green line 
corresponds to the choice j = 2 and £ = 1, ans, finally, red line cor- 
responds to the choice j — £ ^ 1, that is both graphs correspond 
to antisym,m,etrical solutions. 



5.2. Symmetric and antisymmetrical solutions. In order to find symmetric 
and antisymmetrical solutions we set 

q = qi= qi and p = (72 = 93 
In such a case it follows that m = j, that is (p4l) reduces to 



I 2q + 2p = 1 



1,2. 



(35) 



where to j = 2 corresponds sym,m,etric solutions^ while to j = 1 corresponds anti- 
symmetrical solutions. 

We see that the exchange j ^ 3 ~ j reduces to the same system provided that 
1] -^ — ry, E — > —E and £ — >■ 3 — ^; thus we can choose, for argument's sake, j = 2 



14 



ANDREA SACCHETTI 



obtaining 



+ i]p° 



E + U-^Y 



^ = 1,2. 



2q + 2p 



= 1 



where we set 



n = E 



\{-lf 



It immediately follows that the exchange t ^^ i — t reduces to the same system 
provided we perform the exchange p ^ q and E -^ E + ■^{—If. Therefore we can 
choose, for argument's sake, £ = 2 and the system takes the final form 



i-1 






1 



= E- 

= E + 



(36) 



By means of a straightforward calculation we can obtain 77 and Q as functions 
on q and plot D, versus rj. 



(d) 



(c) 



(b) 



(a) 



Figure 5. Here we plot the absolute value of the symmetric wave- 
functions corresponding to the case j = £ = 2 for ?] = — 12 and 
77 = 10, in the case of cubic nonlinearity. For 77 = —12 < —77^ we 
have a wavefunction (a), which is the continuation of the unper- 
turbed one, plus two new wavefunctions (b) and (c) associated to 
branches coming from a saddle point; for 77 = 10 we only have the 
wavefunction (d), which is the continuation of the unperturbed 
one. 
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In Figure S] we consider the value of fl, as function of 77, corresponding to sym- 
metric/antisymmetrical stationary solutions in the cubic case, where a — 1. We 
may see that there exists a critical value 77^ = 8.324 such that for any \ri\ < rf^ we 
only have 4 symmetric/antisymmetrical stationary solutions as in the linear case 
(where 77 = 0). At 7/ = ±7;^ saddle points occur and new branches of symmet- 
ric/antisymmetrical stationary solutions arise. The points denoted by (a), (b), 
and (c) correspond to values of 57 for 77 = — 12 < —77^ and where j = £ = 2; the 
point denoted by (d) corresponds to the unique value of Vi for 7/ = 10 and where 
i = I ~ 2. In particular we have that (see also Figure [5]) 

(a) 77^-12 r2 = -7.021 9 = 0.010 

(6) 77 = -12 r2 = -5.952 9 = 0.478 , . 

(c) 77 = -12 r2 = -4.534 9 = 0.314 ^ ' 

(d) 77 = 10 r^ = 0.990 q = 0.214 

We may remark that in the limit |7;| — >■ -|-cx), that is for large nonlinearity, then 
the wavefunctions (a) and (b) are, respectively, fully localized on the internal and 
external two wells, while the wavefunctions (c) and (d) are equally distributed on 
the four wells. 

Remark 9. We should remark that the same picture occurs even for other values 
of a; for instance in the case of a = 2 the saddle point occurs at rj = rf = —16.648, 
in the case of a — Z the the saddle point occurs at rj = rf^ = —38.775. In generale, 
computing rysigma for higher values of a the following rule appears: rf^^ jrf ^ 2 
for large a . 



5.3. Asymmetrical solutions. Now, we look for the asymmetrical solutions ([34|: 
actually, we have a family of 8 different systems of equations. In fact, we restrict 
our attention to the branches of asymmetrical solutions connected to the symmetric 
ground state. To this end we choose j = ^ = 7ti = 2. The numerical solutions of 
([M)l in the cubic case (i.e. for cr = 1) are plotted in Figure H) more precisely we 
plot the energy fl as function of 77. As we can see a supercritical bifurcation of 
the symmetric stationary solution occurs at 77 « —2.31, the new branch behaves as 
f2 « 77 for large value of I77I, and the four almost-degenerate eigenfunctions are fully 
localized on one single well (see Table [T]). 

We may remark that (see also Figure [71 left panel) that the bifurcation is of 
the same supercritical kind as in the double well model (where bifurcation occurs 
at 77 = —2); furthermore, as in double well model, we observe (see also Figure [3 
right panel) a subcritical bifurcation of the symmetric stationary solution when the 

value of the nonlinearity power is bigger than the critical value ^.threshold — 2 
obtained by |15) . In adjoint to this spontaneous symmetry breaking effect of 
the ground state, which is the most relevant effect, other spontaneous symmetry 
breaking effect of the higher energy stationary states occur, and also new branches, 
associated to saddle points, of asymmetrical stationary states arise (see Figure [6] 
again) . 

5.4. Ground state solution for large nonlinearity. As discussed above, we 
have seen that for |7;| large enough (actually 77 = —12, as computed in Table [T]) 
the four almost-degenerate asymmetrical solutions, associated to the branch which 
bifurcates from the symmetric ground state, are localized on one single well. This 
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a. 



5 







-5- 



-10 



-15 



-20 




-20 -10 10 



r[ 



Figure 6. Here we plot the energy J7 corresponding to symmetric 
(red lines) and asymmetrical (black lines) solutions versus i] in the 
case of cubic nonlinearity (i.e. cr = 1). Here we choice j = i ~ 
TO = 2. The symmetric ground state bifurcates at 77 = —2.31 and 
new asymmetrical solutions occur. The bifurcation point, enclosed 
in the box, will be zoomed in Fig. [71 



result can be proved by means of a simple asymptotic argument as 77 — >■ — 00. More 
precisely, let j — ni — £ — 2 and let us set 



qi 



1 + 77 ^si and n = 7]{l + T] ^r) 



where si = si(?7) and F = T{ri) will be discussed later. From equation (p4l -5) it 
follows that Qk — 0(77^^), k = 2,3,4; more precisely, from equations (|34l-2). (p4l -3) 
and (|Ml-4') immediately follows that 

92 = ??"^S2, 93 = V^'^ss and 174 = 77^*^54 

where Sk '■= Sk(jl) ^-re such that S4 ~ S3 ~ S2, as 77 — > — cx), and where S2 ~ 1, as 
77 — > —00. Since equations ([511 -1) and ([Ml -5) imply that 



si -f S2 ^ and vs2 + asi — F ^ , as 77 —> —00 , 



hence 



Sl 



-1 and F- 1-cr. 
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n 


qi 


92 


93 


94 


-12.000 (ei) 


0.5-10-^ 


0.007 


0.986 


0.007 


-12.000 (ez) 


0.007 


0.986 


0.007 


0.5-10-^ 


-11.999 (eg) 


0.3-10-'' 


0.5-10-^ 


0.007 


0.993 


-11.999 (64) 


0.993 


0.007 


0.5-10-=' 


o.s-io-'' 


-7.021 (a) 


0.010 


0.490 


0.490 


0.010 


-7.009 


0.2-10"^ 


0.011 


0.487 


0.502 


-7.009 


0.502 


0.487 


0.011 


0.2-10-^ 


-5.979 


0.013 


0.451 


0.069 


0.466 


-5.979 


0.466 


0.069 


0.451 


0.013 


-5.952 (b) 


0.478 


0.022 


0.022 


0.478 


-5.376 


0.013 


0.364 


0.240 


0.382 


-5.376 


0.382 


0.240 


0.364 


0.013 


-4.682 


0.347 


0.091 


0244 


0.317 


-4.682 


0.317 


0.244 


0.091 


0.347 


-4.534 (c) 


0.314 


0.186 


0.186 


0.314 



Table 1. Here we collect all the solutions corresponding to the 
value 1] — —12. Wavefunctions labeled with the letter (a), (b) 
and (c) coincides with the symmetrical ones already computed in 
P7|) : the other wavefunctions are asymmetrical wavefunctions. In 
particular the 4 wavefunctions (ei)-(ei) associated to the ground 
state ri ^ — 12 are fully localized on one single well 




-30 -25 -20 -15 -10 -15 1 



Figure 7. In the left picture we zoom the bifurcation point in the 
cubic case, i.e. a — 1. For large value of a the kind of bifurcation 
may change; in the right picture we plot the bifurcation point for 
(7 = 4. Red lines denote the symmetric stationary solutions, black 
lines denote the branch of the asymmetrical stationary solutions. 



This solution corresponds to the ground state, indeed it minimizes the Hamiltonian 
function (pS)) since J2k1k~^^ — IJ^klk]"^ and \e/(3\ is large enough (because we 
are considering the case |?7| >• 1). 
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Since a similar argument may apply when we choose, as starting point, qk = 
1 + rj^^Sk, for fc = 2, 3, 4, then we have proved the following result. 

Theorem 1. There exists a value rj* > such that the symmetrical stationary 
ground state bifurcates at rj — —rf and in the limit of large focusing nonlinearity, 
that is rj ^ — oo, then and the four almost- degenerate asymmetrical solutions, which 
arise at the bifurcation point, are localized on one single well. 

Remark 10. We can extend this result to any number N of wells; indeed the same 
asymptotic argument applies to the system i32\) where we choose 9^ = 9k+i for any 
fc = 1,2,...,N -1. 

Remark 11. By making use of the same arguments in 2 one may prove that the 
resulting almost-degenerate asymmetrical stationary solutions are orbitally stable; 
however we don't dwell here these details. 

6. Conclusion 

Semiclassical methods turn out to be a very powerfull tools in order to reduce 
a NLS to a finite-dimensional Hamiltonian systems. Indeed, by applying such 
techniques jointly with the Lyapunov-Schmidt reduction scheme, we are able to 
describe the ground state solutions as a superposition of vectors localized on single 
wells, with a rigorous estimate of the error [2]. In particular the Hamiltonian 
system (|25|) we obtain it is explicitely written, it can be reduced and it can be 
studied by means of standard numerical tools. 

In more details we consider the case with A^ = 4 wells and we see that the 
spontaneous symmetry breaking effect, already observed in a double well model, 
similarly occurs; in particular we still observe supercritical bifurcation when the 
nonlinearity parameter is less that a threshold value, for value bigger that such a 
threshold value a subcritical bifurcation occurs. 

A remarkable result is that in the case of large enough focusing nonlinearity then 
the symmetric ground state bifurcates and the new N almost-degenerate stable 
solutions are fully localized on one single well. 
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